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Abstract. By use of a generalized Bloch vector construction, we study 
the decoherence of a system composed of two interacting qubits in a noisy 
environment. In particular, we investigate the effects of correlations in the noise 
acting on distinct qubits. Our treatment of the two-qubit system by use of the 
generalized Bloch vector leads to tractable analytic equations for the dynamics of 
the 4-level Bloch vector and allows for the application of geometrical concepts from 
the well known 2-level Bloch sphere. We find that in the presence of correlated 
or anticorrelated noise, the rate of decoherence is very sensitive to the initial two- 
qubit state, as well as to the symmetry of the Hamiltonian. In the absence of 
symmetry in the Hamiltonian, correlations only weakly impact the decoherence 
rate. 



PACS numbers: 03.67.Pp, 03.65.Yz, 74.78.-w, 85.25.Cp 
1. Introduction 

Loss of coherence, or quantum decoherence, remains the most important obstacle to 
overcome in order to build a quantum computer. Decoherence due to noise or to 
entanglement with uncontrollable degrees of freedom is responsible for the decay of 
coherent superpositions of qubit states. The result is irreversible loss of the quantum 
information required for operation of the device. 

Leading candidates for the design of a quantum computer are superconductor- 
based qubits involving Josephson junctions [T], qubits based on nitrogen vacancies in 
diamond [2, 3 , optically trapped ions [1] and quantum devices including semiconductor 
quantum dots [5HH] ■ It is commonly assumed that the decoherence mainly originates 
from low- frequency noise with spectrum of l//-type created by dynamic defects 
tunneling between two metastable states - so-called two- level systems (TLSs). In 
particular, in Josephson charge qubits the TLSs are formed by traps in the amorphous 
material of the insulating barrier of the Josephson junction, and in the substrate 
required to fabricate the qubit. 

While the effect of fluctuating TLSs on isolated qubits were extensively studied, 
see, e.g., references [9rtl5]. or [16] for a review, it is clear that each of the tunneling 
charges present in the insulating substrate might influence several qubits fabricated 
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on the same chip. As a result, the noise acting on different qubits may be correlated, 
which can manifest itself in modifying the decoherence comparing to the case of many 
uncorrelated sources of the same total intensity. Specifically in charge qubits, we 
might have a situation where trapped fluctuators are located in the space between 
the Cooper-pair-boxes belonging to two different qubits that may lead to correlated 
or anti-correlated noise acting on the two-qubit system. 

Quantum oscillations in superconducting two-qubit systems were demonstrated 
already in 2003 [T7] , and have recently been extended to three qubit systems [T51 HH] ■ 
The effect of correlated or partially correlated noise acting on two qubit systems have 
been studied theoretically in [20-- 27 . However, due to the complexity of the Hilbert 
space of coupled qubits, the efforts have primarily resulted in numerical surveying 
of various situations. In addition a general theory for passive error protection in 
decoherence-free subspaces (DFS), subspaces of the multi-qubit Hilbert space which 
are immune to correlated noise, has been developed [2"5H3"2"] or see [33] for a review. 
The work in this paper does not attempt to extend the general work on DFSs, but 
rather provide an efficient way of analyzing the decoherence for initial states that are 
not necessarily in the DFSs. 

In this work, we present a Bloch vector treatment of the two-qubit decoherence 
problem, generalizing the well known formulation for a single qubit. The 15- 
dimensional Bloch sphere of the two-qubit system, and the subspace of available 
positive definite density matrices, has a much richer mathematical structure than the 
corresponding construction for a single qubit (34 - 36 . Still the formalism is similar 
to the single-qubit problem, and we can take advantage of the familiar geometrical 
concepts developed for the treatment of two levels. 

In order to illustrate the 4-level Bloch vector formulation, we consider a model 
problem of two coupled qubits, Qi and Q2, subject to noise from two channels: 
acting on Qi and &(t) acting on Q 2 . We are able to derive simple analytical formulas 
for the decoherence rate of two qubits as a function of the degree of correlation S c 
of the noise sources. Our results clearly demonstrate a strong impact of correlations 
on the decoherence rate that is very different in distinct subspaces of the two-qubit 
Hilbert space. This effect is sensitive to the symmetry of the two-qubit Hamiltonian 
and in the absence of symmetry, correlations in the noise are not important for the 
decoherence rate. 

The rest of the paper is structured as follows. In Sec. [2] we define our model of 
two coupled qubits in a generalized environment, and discuss qualitative aspects of our 
model using a simplified picture. The details of the initial dephasing of our system, 
assuming stationary fluctuators with unknown distribution of initial states, is studied 
in Sec. [3| The long-time decoherence in the weak coupling limit is studied in Sec. [4] 
Finally, our results are discussed in Sec. [5] 

2. Model 

Following |22I) the Hamiltonian describing two coupled qubits, each at zero bias, where 
the qubits are coupled directly by interaction between the basis-forming degrees of 
freedom, can be written as H = H + Hi (t) where 

H = + (i) 

3=1,2 



Bloch-sphere approach to correlated noise in coupled qubits 



3 



2=1,2 

Here £j(t), j £ {1, 2} are random forces acting on the qubits, Q% and Q%, respectively. 
The random forces might origin from the same source, from different sources or a 
combination, resulting in different degree of correlation in the signal. The degree of 
correlation is given by the expression 

s= ((&(*)-&)(&(*)-&)) (3) 



V((ei(*)-^i) 2 X(6W-e 



2) 



2\ 



where the brackets indicate ensemble averaging over initial conditions and realizations 
of the noise process, and & = 

The eigenstates and eigenvalues of the Hamiltonian ([!]) are [22] : 
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where The states |fcZ), /c, I £ {0, 1} are the eigenstates of the operator o~\ ® erf, 



|ifcjl)± = (|ifc) ± |j7))/2, A±=A 1 ±A 2 , 0± = ,/A 
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7± = V(l±A_/n_)/2, r?± = V(l±A+/n+)/2. 

Having defined our model and its parameters, we present our procedure for the study 
of 4-level systems by use of the Bloch- vector construction. 

2.1. Coordinate transformations 

Noise-induced decoherence leads to deviation from coherent oscillations of the states 
due to the Hamiltonian H . In order to study the effect of noise, it is useful to first 
transform our system to the reference frame where the density matrix is stationary 
in the absence of noise. This procedure allows us to separate decoherence from 
rapid quantum oscillations. The transformation is just the standard transform to 
the interaction picture: 

ft = rf(t)p s (t)U(t), V'(t) = U\t)H x (t)U{t) (5) 

where U(t) = e~^ % t^ Hot , and ps(t) is the density matrix of the two-qubit system 



in the Schrodinger picture. The explicit expression for V'(t) is given by Eq. ( A.l I 
in |Appendix A| 

We proceed by specifying the initial state, p' , of our two-qubit system. Next, we 
perform a similarity transform 

p = S~ 1 p'S, V"(t) = S- 1 V'(t)S ) (6) 

where S is the eigenvector matrix of the initial state p' . The motivation behind the 
similarity transform is the following treatment in terms of coherence vectors on the 
generalized Bloch ball. 

For a two-qubit system, the set E4 of positive definite density matrices of trace 
one is a convex subset of the 4-level Bloch sphere in M 15 , contrary to the 2-dimensional 
case where the set of density matrices are equivalent with the set defined by the Bloch 
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sphere. Specifically, the set of density matrices corresponding to pure states is only 
a 6-dimensional subset of the surface of the 4-level Bloch sphere. For more in-depth 
analysis of the geometry of Bloch vectors in the two-qubit system see [37] . The set E4 
can be parametrized by use of the generators Xi of SU4, which are tabulated in |37| . 
Using this parametrization we write 

p = \ + \Y.^{t)\, v"{t) = \Y,ps)\- (7) 

i=l i=l 

The motivation behind the similarity transform, (|6]), is to obtain the special initial 
state 

1 / 

where all the components, mi, of the Bloch vector expansion vanish, except for the 
component mis = — VW^- The factor y3/2 origins from the normalization condition 
for this particular parametrization of the density matrix p, given by 53i=i \ m i\ 2 — 3/2. 
The transform to the coordinate frame given by ^ greatly simplifies the later analysis 
of the equations of motion for p. This is the geometric equivalent to the conventional 
choice of coordinate system for the 2-level Bloch sphere where the initial state is 
parallel to the z-axis. 



Po = S V S=i-y|Ais = 



/ 
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2.2. Equations of motion on the Bloch sphere 

The time evolution of the density matrix in the interaction picture is given by the 
standard von Neumann equation, 

d ^ = -i[V"(t),p(t)]. (9) 

By use of the Bloch vector parametrization of p(t) and V"(t), given by ([7]), we are now 
able to calculate the time evolution of p(t) in a way that simplifies further analysis 
of the decoherence. The time evolution of the coefficients m{t) of Q can be derived 
from the commutator relation 

[Xi,X j ]=2if ijk X k , (10) 

where fijt are the structure factors of SU4 (a table can be found in We obtain 

the following equations for the components: 

rhi(t) = fijkl3 3 (t)m k (t). (11) 

Due to the similarity transform by S, the initial state is always speci fied by 
m; = for i < 15 and mi 5 = —^/3/2. With this choice we can simplify |ri) by 
writing mi(t) — rrii(0) + ct;(i) and 

rhi(t) = fiji5pj(t)[m 15 (0) + a 15 (t)} + fijkPj(t)a k (t) 

for k 7^ 15. To the first order in the parameter £(£) we obtain the approximation 

ati(t) w fijis0 3 it)m 1B (O). (12) 

This approximation is valid for times sufficiently short, such that all components of 
the Bloch vector tangent to the Bloch sphere due to noise are small compared to the 
component along the initial state. 
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For each individual realization of the stochastic process £j(i), the two-qubit 
system evolves as a pure state on the surface of its Bloch sphere. Since the set of 
density matrices corresponding to pure states is only a 6-dimensional subset of the 14- 
dimensional 4-level Bloch sphere, there can only be six nonvanishing components in the 
equations of motion, ( p~2j ) . Averaged over the noise process £j(t), we obtain diffusive 
dynamics of the Bloch vector in this 6-dimensional tangent plane. The probability 
distribution for the Bloch vector will grow in width in each of the six directions, and 
in general the distribution width might grow faster in some directions. 

Proceeding, as long as we consider times sufficiently short so that the uncertainty 
of the Bloch vector due to noise is small relative to the size of the Bloch sphere, we can 
calculate the decay of the component mi 5 (t) = —y/3/2 + ai^(t) along the Ais-axis, 
the by use of the normalization condition 



ai 5 (t) = J ^ 



\ 
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Averaged over initial conditions and realizations of the noise process we obtain 

< a ^)> w \/IE <«?(*)>. (14) 



14 

where the mean square of the components at is defined by 

t t 

<«?(*)> = J j dhdt 2 ai(ti)ai{t2). (15) 
o o 



2.3. Simple picture of decoherence of a two-qubit system by a single TLS 

Before we move on to the calculations it may be instructive to look at a simple 
model that still contains essential features regarding the sensitivity of our two-qubit 
system to external noise. Assume two (coupled) qubits, Q\ and Q2, both coupled 
to an environment containing only a single quantum object, e.g., a TLS, through 
two channels £1 and £2 as illustrated in Fig. [I] a) and described by the interaction 
Hamiltonian 

H I = (^ 1 al+^ 2 z )E = XE. (16) 

Here E is an operator acting on the quantum object. The evolution of the state |</>s(t)) 
of the environmental quantum object is in general dependent on the state of the two- 
qubit system. The decoherence of the two-qubit system is determined by tracing over 
the degrees of freedom of the quantum object, i.e., by the overlap matrix elements 
(^sWI^eW)- Here i,j denote the four orthogonal states of the two-qubit system. 
In other words, the decoherence is determined by the difference in the conditioned 
evolution of the state \4> l E (t)) depending on the initial state of the two-qubit system. 

In the configuration described in figure [I] b) , charge fluctuations between the 
impurity and the metallic gate give origin to a correlated noise on the two qubits. If 
each qubit couple with the same strength to the gate, £1 = £2, and the two qubits 
are prepared in a superposition of the states 1 01) and 1 10) then the expectation value 
of operator X coupling to the environment will vanish. The environmental quantum 
object is therefore unable to distinguish between these two states and therefore no 
entanglement build-up can take place between the environment and the qubits. In 
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Figure 1. a) Two qubits, Qi and Q2, coupled to a TLS in the substrate between 
the two qubits, where the charge can tunnel between two sites, b) The qubits 
are coupled to a charged impurity through its image charge on the metallic gate, 
the charge can tunnel between the gate and the impurity. The coupling strengths 
are given by §1 and £2- The configuration demonstrated in a) give origin to 
anticorrelated noise, while configuration b) give origin to correlated noise. 

this model, the states 1 01} and 1 10) span a decoherence-free subspace (DFS) in the 
presence of correlated noise [213 [33] . However, if the coupling parameters £1 and 
£2 are not exactly the same the states 1 01} and 1 10} will induce different potentials 
on the environment such that entanglement can build up. In the extreme case, where 
£1 = — £ 2 , which could in fact be realized by a charge trap placed between two charge 
qubits, as could be realized in a configuration similar to figure [l] a), completely 
correlated noise will be twice as efficient as the same noise from two uncorrelated 
sources. Thus this possible situation needs to be accounted for in the design process 
of many qubit cirquits. 

Furthermore, if the two-qubit system has dynamics given by an intrinsic 
Hamiltonian H , initial states that are not eigenstate of H will in general evolve in 
time. In this case, we will have decoherence even if the noise is completely correlated, 
£1 = £,2, and the initial state is 1 01) due to the fact that the state will decompose in 
eigenstates which in general might not be protected from correlated noise. 

3. Initial decoherence in the stationary path approximation 

In this section, we will study the initial decoherence of our two-qubit system focusing 
on the sensitivity of the decoherence to correlations in the noise. Previously it has 
been found that the sensitivity to noise in two-qubit systems is reduced if the noise 
is correlated, and the two-qubit system is prepared to be in a state within the SWAP 
subspace, spanned by the states {|10) , |01)} [22] . In order to check this, we calculate 
the decoherence rate from four different initial states, two of which belongs to the 
SWAP subspace, and two belonging to the subspace spanned by 1 00) and |11). 

As an example to illustrate the method we will calculate the initial decoherence 
rate when the two-qubit system is initially prepared in the singlet state |^o) = 
1 10) — 1 01), where we have omitted the normalization factor. Proceeding as described 
in section [2] we transform to the coordinate frame where the initial Bloch vector of 
the two-qubit system lies on the south pole of the Bloch sphere, m, = —a/3/2 for 
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15 and otherwise. The transformation matrix is for this initial state given by 



S 



1 



/ V2 



^ 
0-10-1 
0-101 
V V2 J 



(17) 



We can then compute the noise matrix 



1 



15 



V" =S T V'S= 

(=1 



and use (12) to compute the coefficients a(t) parameterizing the density matrix p(t). 
For the initial state \?po} = 1 10) — 1 01) and the Hamiltonian given by ([2|, the differential 
equations for en are given by Eq. (A.2). We note that there are only six nonvanishing 
components on. 



Since the general formulas given by (A.2) are rather complicated, it is instructive 



to investigate the simplified case where the two-qubit system is at the co-resonant 
point Ai = A2 = A. The problem simplifies further if we assume that the inter-qubit 
coupling is much less than the single qubit tunneling element, v <C A. With this 
symmetric Hamiltonian, we get a simplified set of parameters j± = 1/^/2, rj± = 0. 

Furthermore, if we are interested only in the initial decoherence rate, we may 
assume the stationary path approximation, where the noise sources are stationary 
and decoherence is solely due to the different realizations of initial values of the noise 
sources at the start of each experiment (e.g., a statistical distribution of the gate 
charge at the start of each run of the measurement protocol [25 ). After averaging 
over many experiments with different initial conditions for the noise sources we obtain 
an uncertainty in the precise value of the Bloch vector of our two-qubit system. By 



use of (15), inserting the expressions for the derivatives of the component Ofj given in 



(A.2) we obtain the following expressions for the mean square deviations along each 



component due to noise: 

(«i(*)> 
<«§(*)> 



(c4(*)> = (1/8) ((j£) 2 _> m*) ■ 

(a l3 (t)) = am(mt)is-(t) 



-(*)] 



where (5£)z 



(a 2 n (t)) =(l/4)((^) 2 _)[c_(i)-c + (t)] 2 : 
(<&(<)) = (1/4) ((SO 2 -) [s-(t) + s + (t)} 2 

6(0) t 6(o), w± = n + ±n_, s± (t) = 



L (i 



(18) 

cosw±t), s±(t) = 



u>± smuj±t. The expression is valid for times shorter than the correlation time of the 
noise sources, t <C t c . 



From the normalization condition, (14), and the expressions for the mean square 



fluctuations along the six components, (18), we obtain the following simple formula 



determining the initial decoherence of the two-qubit system: 



(«16(*)) 



<(^) 2 -> 



/(*), /(*) 



C_(i) 



-(*) 



(19) 



V6 w_ lu + 

As expected, we find that at the symmetric co-resonance point the coherence of the 
initial singlet state 1 10) — 1 01) is very sensitive to the degree of correlation in the 
noise sources 6 and £2- If t ne noise sources are completely correlated there is no 
decoherence. 
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(a 15 {t)} ■ 


|01) - |10) 
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|01) + |10) 


(m 2 -)[f(t)+9(t)] + ((S0l)[f(t)-9(t)} 


|00) + |11) 




|00) - |11) 





The expressions for (ai§(t)) for different initial states are summarized in the 
Table [1] In this table, 

g(t) = (w + o;_) _1 (l — cosw_i — cosu + t + cos2f2i). (20) 

In the absence of correlation 

mt) = mi) = 2 ((Kiwi 2 ) + (mo)} 2 )) , 

and for the fully correlated noise ((5£)?_) — 0, while for the anti-correlated noise 

m)i)=o. 

We find that the effect of correlated noise is strongly dependent on the initial 
state of the qubit. The singlet state \tp) = 1 01) — 1 10) is a decoherence-free state in 
the presence of completely correlated noise, while the state 1 01) + 1 10) is only partly 
protected from correlated noise. The difference is due to the fact that the singlet state 
is an eigenstate of Hq while the the initial state 1 01) + 1 10) is not an eigenstate of 
Hq and will therefore due to the time evolution obtain components in subspaces of 
the Hilbert space that are sensitive to correlated noise. For a general initial state, 
|(/>o), the component in the subspace of the Hilbert space spanned by the singlet state, 
( (01 1 — (10|)|</>o), will not decay in time. For this state we will find an initial decay 
of the coherence, and then persisting coherence of amplitude given by the overlap 
element. The initial states \ipo) = 1 00) — |11) and |-0o) = 1 00) + |11) are decoherence- 
free and weakly protected, respectively, to anticorrelated noise. It is important to 
notice that we have assumed co-resonance. If the Hamiltonian, ([!]), is not perfectly 
symmetric the symmetric states will no longer be eigenstates of Hq and will therefore 
be less sensitive to correlated/anticorrelated noise. 

To compare with the above result, if the two-qubit Hamiltonian commutes with 
the operator representing coupling to the noise source, i.e., for the conventional 
Hamiltonian [35] 

J=l,2 3=1,2 

then the subspace spanned by the set of states 1 01) , 1 10) is a 2-dimensional DFS with 
respect to a completely correlated noise, but vulnerable to an anticorrelated noise, 
and vice versa for the set of states spanned by |00) , |11). 

4. Decoherence for intermediate times in the weak coupling limit 

In the previous section we studied the initial dynamics of the two-qubit system, due to 
noise from the channels and £2(^)5 by use of the stationary path approximation. 



Bloch-sphere approach to correlated noise in coupled qubits 



9 



In this section we will consider the dccohercncc for times longer than the correlation 
time t c of the noise sources. 

We proceed by use of the four-level Bloch-vector method applied in the same 



manner as presented earlier, by (111, (14) and (151. The decoherence is restricted to a 
six dimensional subspace of the 14-dimensional Bloch-sphere as long as we observe the 
qubits for times 1 1£| <C 1. We are required to solve the integrals, of the form given by 



(15 1, determining the mean square component in each of the six directions spanning 
the subspace. 

To compare with the results of the references [HIES], we start with the initial state 
\ipo) = 1 01) . After transforming to the coordinate frame where p — 1/4 — \/3/8Ai5, 
we obtain the following equations of motion 

&9(t) = £i(t)(\i sin cu_t + \ 2 sin uj + t), 
&io(t) = £i(*) (/•*!- + A*2-)(cosw_i - cosu;_i), 
au(t) = [fi(t) - 6(f)] [-(M2+ + sina;_t 

+(/i 2 + - fi 1+ )sinuj + t] , 
d 12 (i)=0, (21) 
"13(f) = ?2(0(-A 2 sinw_t - Ai sino; + f), 
«14(f) = £2(^X^1- — A*2-)( — cosw_t + COSCJ+f) 

where we have introduced the notation 

\ 22 22 \ 22 22 

Ai =ry_7_ -7j + 7 + , A 2 = - ry + 7_, 

Mi± = »7+»?- (7+ ±7-)) A^2± = 7+7-0?- ± «7+)- 

When computing the mean square of each component (af(£)), by use of (15), 
we assume time translation invariance, which allows us to make the standard 
transformation 140 



(«?(*)) = J dT J dra t (- T - + T) d, g + Tj , (22) 
-00 

where r = t 2 — ti and T = (ti + f 2 )/2. We have also assumed observation times to be 
longer than the correlation time of the noise sources, t 3> t c , which allows us to extend 
the integral over r to infinity. In addition, we have to assume that the coupling to 
the environment is sufficiently weak to ensure that the six components of the initial 
"hypersurface" where dynamics take place, (af(t)) for i = 9 — 14, are small compared 
to (m 15 (t)) . 



Inserting (22) into (22), we obtain terms in the integrand proportional to 
cos(J7 + ± f2-)r, cos 2(0+ ± 0_)T and cos(±il±r ± 2f2±T). If we assume that 
i > 1/ cj + , the terms that oscillate as a function of T will be negligible compared to 
the slowly varying ones. By this approximation the expressions for the mean square 
of the six components, from the initial state |^o) = 1 01) , reduce to the following: 



t 00 

)(*)) rfr (^( r )^(°))[ A ? COSW - r + A 2 COSW + T ]- 



(23) 



Here we omitted the explicit expressions for the other five components. We notice 
that the integral grows linearly in time. 
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By use of the normalization condition, (14 1, we can make the following 
approximation for t |£i j2 | *C 1, 

(a 15 (t)) « (1/V6) $>?(<)) « ff/V6 (24) 

where we introduced the decoherence rate T characterizing the decay of the Bloch 
vector given by (a^i)) to lowest order in t. With the initial condition \ipo) = 1 1 ) we 
obtain the following expression for T: 

r = (i/Vq){ Y, SuK) + s 22 {u ± )] a± 
± 

- 2S 12 {uj-){h2+ + Mi+) 2 - 2S 12 (uj + )(jp 2+ - A*i+) 2 }, 
a + = \\ + {p 2+ + fi 1+ ) 2 + (ju 2 - + Mi-) 2 , 

a- = A 2 , + ([i 2+ - Mi+) 2 + (Ma- + Mi-) 2 ■ (25) 
Here we have introduced the spectral density 



Sy(u;)= / dr(6(r)0(0))cos( W T) (26) 

— oo 

for the Fourier transform of the (cross) correlation spectra from channels i,j at 



frequency uj. The expression given by (25) is valid for max{l/w + , r c } <C t <C 1/ |£i,2| 



We note that the noise at the two frequencies, f2+±f2_, is the only relevant component 
of the noise in the weak coupling limit. This result agrees with the findings of [22]. It 
is, however, important to stress that the above result is to first order in the coupling 
strength For a single qubit, with energy splitting A, subject to transverse noise, 
it is known [41] [42] that we will get a contribution to the dephasing from the noise at 
zero frequency 

(a^t)) cx£ 2 SH+£ 4 S(0) + .., 

where uj — A is the qubit frequency at zero bias, and the zero frequency contribution 
is due to the net shift in the precession rate due to the noise, uj = ^/A 2 + £ 2 = 
A + £ 2 /2A 2 . This effect will also be present in the multi-qubit problem. 

In the work by Averin and Rabenstein [22], it was also reported that the 
decoherence rate was more sensitive to correlation when the two-qubit system was 
prepared in the state \ifjo) = 1 01) than in the state \ipo) = 1 00) . Contrary to these 
results, we find that the two states are initially equally sensitive to correlations. If 
the noise sources have the same magnitude, Sii(u;) = S 22 (oj), we can express the 
decoherence rate as a function of the correlation parameter Sc given by ([3]). The rate 
as a function of correlation is plotted in figure [2j where we find that the coherence 
time of the state jV'o) = 1 1) is enhanced in the presence of correlated noise, while the 
coherence time of the initial state |^>o) = 1 00) is diminished by the same factor. Our 
results is fully consistent with those of Averin, and the "seeming" discrepancy has the 
following explanation. In reference [22], the oscillations of the probability pi, 2 (t) to find 
Qi and Q 2 in the excited state at time t are plotted. Our results for the decoherence 
rate agree at short times (in the regime t\£\ <C 1 where our approximation is valid), 
but thereafter in [22] slowly damped oscillations of amplitude ~ 0.7 for the initial state 
1 01) were found. This behavior is due to the fact that the state 1 01) is not an eigenstate 
of Hq. For the Hamiltonian (JlJ at the co-resonance point, the eigenstates are the usual 
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A 1= 1.0 




s s 

c c 

Figure 2. (Color online) Sensitivity of the decoherence rate T to correlations 
between £i(t) and $2(t) and to symmetry in the Hamiltonian Ho. We plot the 
total decoherence rate V against the correlation coefficient Sc for four different 
initial states, |V>o} = [01) - 1 10) (red), \tp ) = |01) (blue), \ip ) = |00) (black) and 
\ipo) = 1 00) — 1 1 1) (green). In the top figure we are at the co-resonance point, 
where Ai = A2, while in the bottom figure Ai = O.7A2. The coupling parameter 
is v = 0.1. 

singlet and triplet states. The singlet state is protected against correlated noise, while 
the triplet is not. The initial state 1 01) has equal components in both subspaces. In 
the presence of a correlated noise, the component in the triplet will therefore decay 
rapidly while we will find persisting oscillations due to the finite component in the 
singlet subspace. The Hamiltonian considered in |22j is not fully symmetric, but the 
picture is still valid. One component decays at a fast rate, while the component in 
the approximate "singlet" -subspace decays much slower. Interestingly, it was reported 
in that concurrence oscillations from the initial state 1 01) lasts even longer than 
expected from its components in the DFS and its compliment, due to an interference 
effect between the components that decreases the decoherence. 

From a practical point of view, however, we can focus on the initial decoherence 
rate. If the purity of the state has decayed more than a few percent quantum 
computation by error correction is practically impossible. Furthermore, our results 
show that the initial state is even more important than it was previously reported 
[22"1 125] , In the presence of completely correlated noise, Sc = 1, the singlet state 
\4>o) = 1 01) — 1 10) is decoherence-free. However, it is important to mention that 
this is the only decoherence-free state. There is no 2-dimensional DFS for our 
Hamiltonian Hq. The state \tpo) = 1 01) — 1 10) is decoherence-free if and only if 
we are at the co-resonance point, as shown in figure [2j If we move away from this 
point the sensitivity to correlations is reduced. This behavior is controlled by the 
symmetry of the Hamiltonian. For our Hamiltonian, ([!]), the relevant parameters are 
Ki = A_/A + and k 2 = A + /^. For either Kx -C 1 or k 2 <C 1 the Hamiltonian has a 
high degree of symmetry and there exist states that are protected against correlated 
or anticorrelated noise. If the Hamiltonian has low degree of symmetry, initial states 
that are initially symmetric or antisymmetric with respect to the noise source, i.e., 
the state \ipo) = 1 01) — |10), will decompose in the non-symmetric eigenstates of H 
and quickly lose its symmetry by the time evolution induced by H . In this case the 
sensitivity to correlations in the noise source will be very weak. 
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5. Discussion 

Our Bloch-vector treatment of the two-qubit decoherence problems shows that the 
two qubit problem can be separated in two regimes. In the symmetric regime, the 
symmetry of the initial state |^o) is conserved by the Hamiltonian in the absence of 
noise. In this regime, the system is very sensitive to correlations or anticorrelations 
in the noise sources £i(t) and £a(*) acting on qubits Qi and Q2, respectively. This 
sensitivity, which is consistent with the general DFS picture, is dependent on the 
initial state and on the parameters of the problem in a more complex way than it 
was reported in [2H [53] , where it was assumed that states in the SWAP subspace 
spanned by the states 1 01) and 1 10) should be less influenced by a correlated noise. At 
the same time, the decoherence of states outside this subspace should only be weakly 
influenced by correlations. This thought is based on the concept of DFSs defined as 
the states where the the dissipative part of the Markovian master equation is zero |32) , 
but does not take fully into account the intrinsic dynamics of the system. A general 
theory for dynamically stable DFSs was developed in [32] • We find that it is both the 
symmetry of the initial state and how much this state overlaps with an eigenstate of 
the Hamiltonian in the absence of noise that determines the rate of decoherence. For 
the Hamiltonian given by ([2|, at the co-resonance point, all the three initial states 
1 01) , 1 10) and 1 01) — 1 10) look similar initially for an external noise source coupled 
diagonally in the basis-forming states, see figure [l] However, only the last state is an 
eigenstate of H and therefore fully protected during the time evolution, i.e., a true 
DFS. 

While previous works, except for the general classifications of DFSs, have mainly 
focused on a correlated noise in the SWAP subspace, we find that this is not the only 
interesting subspace in this regard. Considering two qubits at the co-rcsonance point 
we find that the subspace spanned by the states 1 00) , 1 1 1) is protected against an 
antisymmetric noise. For this system, an antisymmetric noise also leads to increase of 
the initial decoherence rate T for states in the SWAP subspace, and an increase in T 
by factor 2 for the singlet state. 

For the case of noise due to charged traps in the amorphous substrate, a trap 
located between two qubits, giving rise to anticorrelated noise, seems just as plausible 
as the opposite situation of strict correlations. In an experimental setup, one might 
take advantage of this situation by first looking for correlated or anticorrelated noise, 
and thereafter decide the encoding of the two qubits in the subspace that is less 
sensitive to the noise present in the relevant setup. 

Our treatment, by use of the 4- level Bloch vector construction, gives a tractable 
set of equations compared to alternative methods. The familiar geometrical concepts 
used when treating two-level systems by the same method can be carried over to 
higher dimensions. In our treatment, different rates of decoherence along different 
directions on the Bloch sphere, as well as their sensitivity to correlations in the 
noise, follow naturally from the formalism. One may hope that this method might 
be extended also to even higher level systems, e.g., to multi-qubit systems relevant 
for the analysis of the noise in a quantum computer. While we are able to derive 
simple analytical expressions for the two-qubit system, the three-qubit system will 
require a 24 dimensional Bloch-sphere, and obtaining analytical expressions in this 
case seems cumbersome except in very special limits, due to the diagonalization of the 
system Hamiltonian Hq and of p, which are both 8x8 matrices. The method is still 
suitable for numerical impementation, where one might take advantage of the familiar 
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geometric concepts from the two-level problem. 

In conclusion, we have demonstrated a generalized Bloch-sphere method by 
treating the problem of two coupled qubits coupled to statistically-independent or 
correlated noise sources. By use of this method we find expressions for the dephasing 
of the two-qubit system as a function of the degree of correlation of the noise sources. 
At the co-resonance point, the decoherence rate is very sensitive to the initial two- 
qubit state. In the absence of symmetry of the Hamiltonian, correlation in the noise 
seems to have very weak impact on the decoherence rate. 

Appendix A. Detailed expressions 
Appendix A.l. The noise matrix V (t) 

The explicit expression for the matrix elements of V'(t) = is given by: 

Ki = £i[(f?-7- + V+J+) 2 cosw_t + (??_7 + - ?]+ 7_) 2 cosw+t] 

+ 6[(f7-7+ +??+7-) 2 cosw_t + (?7_7_ - 77 + 7 + ) 2 cosw+i], 



V L = (6 +6)H(M2+ +^i+)sinu;_t + i(p-2+ - Mi+) sin 
V22 = £i[0?-7- + V+1+) 2 cosu-t + (77.7+ - ?7 + 7_) 2 cosw+t] 

-&[(»7-7+ + cos^-t + (V-1- - V+l+f cosw+t], 

v 23 = (6 _ (P2+ + Mi+) sinw_t - i (/z 2 + - Mi+) sinw+i], 
^24 — £i[ — coscj_t + /U2-) + coscj + t + /t*2-) — isinu;_/;Ai — i smu + tX2], 
V 33 = - £i[(V-l- + ri+l+ f cosw_i + (r/_7 + - 77+ 7 _) 2 cosw+t] 

+ 6[(j7-7+ +?7 + 7-) 2 cosw_t + (r/_7_ - ry + 7 + ) 2 cos u+t], 
V34 = ^[(Mi- — M2-) cosu;_t — (/ii- — M2-) coscj + t — «A 2 sinaj_rj — iX\ sinw+i], 
V44 = -£i[(»7-7- +?7+7+) 2 cosw_f + (?7_7 + - r? + 7_) 2 cos w+i] 

-&[(»7-7+ +?? + 7-) 2 cosw_f + (ry_7_ - ?7 + 7 + ) 2 cosw + t], (A.l) 



where we have only given the upper triangle of the Hermitian matrix V. 

Appendix A. 2. Differential equations for the density matrix components on 

The differential equations for the components of the density matrix evolving from the 
initial state |V>o) = |10> — |01): 




^2[( — Mi- + M2-) coso;_f + (fii- — H2-) cos u + t — i\ 2 smuj-t — iXi sincj_t] 
£i[(Mi- + M2-) cosu;_i — (fi\- + (12-) cosuj + t — iX\ sincj_t — i\ 2 sinw+t], 



aw 



[Ai sinw_i + A 2 sin u; + i] — — j= [X 2 sinu;_t + Ai sin 
[ - + M2-) cosu;_i + + /i 2 _) coscj+t] 



aio 



V2 



V2 




"12 



COS w + t] 
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+ 



[('7-7+ + 1 1+1-) 2 cosuj_t + (?7_7_ - ??+7-) 2 cos 



OJ + t], 



V2 



<*13 = 



aw 



[Ai sinw_t + A 2 sin u; + i] — [A 2 sinu;_t + Ai sin 

v2 

[ — (^1- + /i2-) cosu;_i + + /i 2 _) cosaj + t] 



V2 
6W 



«14 = 



V2 



+ 



[ - (Ail- - M2-) COSLU-t + - ,U 2 _) COS 



u+t]. (A.2) 
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